GENERALIZED WILCZYNSKI INVARIANTS FOR 
NON-LINEAR ORDINARY DIFFERENTIAL 
EQUATIONS 

BORIS DOUBROV 

Abstract. We show that classical Wilczynski-Se-ashi invariants 
of linear systems of ordinary differential equations are generalized 
in a natural way to contact invariants of non-linear ODEs. We 
explore geometric structures associated with equations that have 
vanishing generalized Wilczynski invariants and establish relation- 
ship of such equations with deformation theory of rational curves 
on complex algebraic surfaces. 



1. Introduction 

This paper is devoted to a very important class of contact invariants 
of (non-linear) systems of ordinary differential equations. They can be 
considered as a direct generalization of classical Wilczynski invariants 
of linear differential equations, which are, in its turn, closely related to 
projective invariants of non-parametrized curves. 

The construction of these invariants is based on the fundamental 
idea of approximating the non-linear objects by linear. In case of dif- 
ferential equations we can consider linearization of a given non-linear 
equation along each solution. Roughly speaking, this linearization de- 
scribes the tangent space to the solution set of the given equation. Note 
that in general this set is not Hausdorff, but we can still speak about 
linearization of the equation without any loss of generality. 

Unlike the class of all non-linear equations of a fixed order, which is 
stable with respect contact transformations, the class of linear equa- 
tions forms a category with a much smaller set of morphisms. This 
makes it possible to describe the set of all invariants of the linear equa- 
tions explicitly. The main goal of this paper is to show the general 
geometric procedure that extends these invariants to the class of non- 
linear equations via the notion of the linearization. 
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In short, the major result of this paper is that invariants of the 
linearization of a given non-linear equation are contact invariants of 
this equation. 

In fact, these invariants have been known for ordinary differential 
equations of low order: third order equations [3, 14], fourth order 
ODEs [1, 8, 6], systems of second order [7, 9], but it was not under- 
stood up to now that these invariants come from the linearization of 
non-linear equations. 

In the theory of linear differential equations Wilczynski invariants 
generate all invariants and, in particular, are responsible for the trivi- 
alization of the equation. In non-linear case the generalized Wilczynski 
invariants form only part of the generators in the algebra of all contact 
invariants. However, the equations with vanishing generalized Wilczyn- 
ski invariants have a remarkable property, that its solution space carries 
a natural geometric structure (see [1] and [6]). We discuss this struc- 
ture in more detail in Section 5 for a single ODE. Finally, in Section 6 
we show that all results of this paper can be extended to the case of 
systems of ordinary differential equations. 
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of partial differential equations". I would like to thank organizers of 
this workshop for inviting me to give these lectures. I also would like to 
thank Eugene Ferapontov, Rod Gover and Igor Zelenko for corrections 
and valuable discussions on the topic of this paper. 



As an example, let us outline this idea for the case of a single non- 
linear ODE of order > 3. Note that all equations of smaller order are 
contact equivalent to each other and, thus, do not have any non-trivial 
invariants. 

The complete set of invariants for a single linear ODE was described 
in the classical work of Wilczynski [18]. Namely, consider the class of 
linear homogeneous differential equations 



In fact, these are the most general transformations preserving the class 
of linear equations. 

A function / of the coefficients Pi(x) and their derivatives is called 
a relative invariant of weight I if it is transformed by the rule / i— > 



2. Naive approach 
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(A')'J under the change of variables (2). In particular, if such relative 
invariant vanishes identically for the initial equation, it will also vanish 
for the transformed one. Relative invariants of weight are called 
(absolute) invariants of the linear equation. 

E. Wilczynski [18] gave the complete description of all relative (and, 
thus, absolute) invariants of linear ODEs of any order. It is well-known 
that each equation (1) can be brought by the transformation (2) to the 
so-called Laguerre-Forsyth canonical form: 

y {n+1) + q n -2{x)y {n - 2) + ■■■ + q (x)y = 0. 

The set of transformations preserving this canonical form is already a 
finite-dimensional Lie group: 

/ ax + b ey \ 

This group acts on coefficients q Q , . . . , q n _ 2 of the canonical form, and 
the relative invariants of this action are identified with the relative 
invariants of the general linear equation. The simplest n — 1 relative 
invariants 3 , . . . , 9 n+ i linear in q\^ have the form: 
(3) 

^-2^ L ) (k-j)\(j-i)\ n ~ k+j «-a,...,n + i. 

Each invariant 9 k has weight k. Wilczynski proved that all other invari- 
ants can be expressed in terms of the invariants 9k and their derivatives. 

Although these formulas express invariants in terms of coefficients 
of the canonical form, they can be written explicitly in terms of the 
initial coefficients of the general equation. Moreover, it can be shown 
that each of these invariants 0j is polynomial in terms of functions Pi(x) 
and their derivatives. 

Equivalence theory of linear ODEs is intimately related with 
the projective theory of non-parametrized curves. Namely, let 
{yo(x), . . . ,y n (x)} be a fundamental set of solutions of a linear equa- 
tion S given by (1). Consider the curve 

Le = {[yo(x) : yi(x) : • • • : y n (x)] | x E JR} 

in the n-dimensional projective space MP™. Since the solutions do not 
vanish simultaneously for any this curve is well-defined. More- 

over, since the fundamental set of solutions is linearly independent, 
this curve is non- degenerate, i.e., it is not contained in any hyperplane. 



4 



BORIS DOUBROV 



Finally, since the set of fundamental solutions is defined up to any non- 
degenerate linear transformation, we see that the curve L £ is defined 
up to projective transformations. 

Consider now what happens with this curve, if we apply the transfor- 
mations (2) to the equation £. The transformations (x, y) \— > (x, fJ>(x)y) 
do not change the curve, since they just multiply each solution by n(x). 
The transformations (x, y) \— > (X(x),y) are equivalent to reparametriza- 
tions of L £ . 

Thus, we see that to each linear equation £ we can assign the set 
of projectively-equivalent non-degenerate curves in WP n . It is easy 
to see that this correspondence is one-to-one. Indeed, having a non- 
degenerate curve in MP n we can fix a parameter x on it and write it 
explicitly as [yo(x) : yi(x) : ••• : y n (x)], where the coordinates yi{x) 
are defined modulo a non-zero multiplier fi(x). Since the curve is non- 
degenerate, the set of functions {2/j(x)}j=o,...,n is linearly independent 
and defines a unique linear equation £ having these functions as the 
set of fundamental solutions. 

Each relative invariant / of weight k can be naturally interpreted as 
a section of the line bundle S k (TL) invariant with respect to projective 
transformations. In particular, relative invariants of weight are just 
projective differential invariant of non-parametrized curves in the pro- 
jective space. First examples of such invariants were constructed by 
Sophus Lie [12] for curves on the projective plane and then generalized 
by Halphen [10] to the case of projective spaces of higher dimensions. 
Note that they can be constructed via the standard Cartan moving 
frame method. 

E. Wilczynski also proved the following result characterizing the 
equations with vanishing invariants: 

Theorem 1. Let £ be a linear ODE given by (1). The following con- 
ditions are equivalent: 

(1) invariants 9 3 , . . . , 9 n+1 vanish identically. 

(2) the equation £ is equivalent to the trivial equation = 0. 

(3) the curve Lg is an open part of the normal rational curve in 
RP n . 

(4) the symmetry algebra of L £ is isomorphic to the subalgebra 
s[(2,R) C sl(n + 1,K) acting irreducibly on R n+l . 

Let us show how to extend Wilczynski invariants to arbitrary non- 
linear ordinary differential equations via the notion of linearization. 
Indeed, consider now an arbitrary non-linear ODE solved with respect 
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to the highest derivative: 

(4) y^ = f(x,y,y',...,y^). 

Let Uq{x) be any solution of this equation. Then we can consider the 
linearization of (4) along this solution. It is a linear equation 

where all coefficients are evaluated at the solution yo(x). It describes 
all deformations y € (x) = yo(x) + eh(x) of the solution yo(x), which 
satisfy the equation (4) modulo o(e). 

Consider now Wilczynski invariants 9 3 , . . . , 9 n+ i of the lineariza- 
tion (5). They are polynomial in terms of the coefficients and 
their derivatives (evaluated at the solution y (x)). In general, let 
F(x, y, y\ . . . , y^) be any function of x, y(x) and its derivatives eval- 
uated at any solution of the equation (4). Then, differentiating it by x 
is equivalent to applying the operator of total derivative: 



dx dy ch/™ _1 ) dy^ 

Thus, analytically, Wilczynski invariants of (5) can be expressed as 
polynomials in terms of {^r^ (0 < i < n, j > 0) and are indepen- 
dent of the solution yo(x) we started with. More precisely, there are 
(n — 1) well-defined expressions W 3 , . . . , W n +i polynomial in terms of 

Di {jf^j^Ji which, being evaluated at any solution yo(x) of the equa- 
tion (4), give Wilczynski invariants of its linearization along yo(x). We 
call these expressions W 3 , . . . , W n+ i the generalized Wilczynski invari- 
ants of non-linear ordinary differential equations. 

To understand, what geometric objects correspond to generalized 
Wilczynski invariants we need to consider the jet interpretation of dif- 
ferential equations. Let J % = J' l {E?) be the space of all jets of order i 
of (non-parametrized) curves on the plane. Then for each curve L on 
the plane we can define its lift to the jet space J 1 of order i. Then 
the equation (4) can be considered as a submanifold £ C J n+1 of codi- 
mension 1, and, by the existence and uniqueness theorem for ODEs, 
all lifts of its solutions form a one-dimensional foliation on S. Let us 
denote its tangent one-dimensional distribution by E. Let (x, y) be any 
local coordinate system on the plane. Then it naturally defines a local 
coordinate system (x, y , yi, . . . , yi) on J 1 such that the lift of the graph 
(x,y(x)) is equal to (x, y(x), y'(x), . . . , y®). In these coordinates the 



6 



BORIS DOUBROV 



submanifold £ is given by equation y n +i = f(x, yo, yi, ■ ■ • , y n ), func- 
tions (x,y ,yi, . . . ,y n ) form a local coordinate system on £ and the 
distribution E is generated by the vector field 

„ d d d d 

dx dy dy n „i dy n 

The expressions W3, . . . , W n+ \ can be naturally interpreted as functions 
on the equation manifold £. In fact, each invariant Wi, being a relative 
invariant of weight % of all linearizations, defines a section Wj of the 
line bundle S l E* by the formula: 

Wi(D, D,...,D)=Wt. 

The largest set of invertible transformations preserving the class of 
ODEs of fixed order % consists of so-called contact transformations, 
which are the most general transformations of J 1 preserving the class 
of lifts of plane curves. The main result of this paper can be formulated 
as follows: 

Theorem 2. The sections Wi, % — 3, . . . ,n + l of the line bundles S l E* 
are invariant with respect to contact transformations of J n+1 . 



3. Algebraic model of Wilczynski invariants 

In this section we give an alternative algebraic description of 
Wilczynski invariants, which is due to Se-ashi [15, 16]. There are two 
main reasons for providing an algebraic picture behind Wilczynski in- 
variants. First, it can be easily generalized to invariants of systems of 
ODEs or even to more general classes of linear finite type equations 
(see [15]). Wilczynski himself described only invariants for a single 
ODE of arbitrary order and for linear systems of second order ODEs. 
It seems that analytic methods become too elaborate to proceed with 
systems of higher order. 

Second, Se-ashi's construction gives alternative analytic formulas for 
computing Wilczynski invariants, which are independent of Laguerre- 
Forsyth canonical form. In particular, this explains why Wilczynski 
invariants are polynomial in the initial coefficients and their derivatives, 
while the coefficients of the canonical form are not. 

Denote by V n the set S n (R 2 ) of all homo geneous polynomials of de- 
gree n in two variables i>i,i>2- The standard GL(2, IR)-action on M 2 
is naturally extended to V n and turns it into an irreducible GL(2,M)- 
module. Denote by p n : GL(2,M.) — > GL(V n ) the corresponding repre- 
sentation mapping. 



GENERALIZED WILCZYNSKI INVARIANTS FOR NON-LINEAR ODES 7 



Consider the corresponding action of g((2,R) on V n = S n {R 2 ). De- 
note by X, Y, H, Z the following basis in $1(2, R): 

x = (Vo),H = a» 1 ),Y = (r ),z = (} ) °). 

Then the action of these basis elements on V n is equivalent to the action 
of the following vector fields on R^ff!,^]: 

9 d d TA d d d 

X = V 1 — , H = V l - V 2 TT- , y=V2TT- , Z = V ± - h V 2 — ■ 

OV2 OV\ OV2 OVi OVi OV2 

In the sequel we shall identify $j[(2,R) with its image in gl(V n ) defined 
by this action. 

We define a gradation on V n such that the polynomials E — u™, E\ — 
f"~ 1 f2, . . . , E n — v 2 have degrees — n — 1, —n, . . . , —1 respectively. 
Then the elements X, H, Y, Z define operators of degrees —1, 0, 1 
and respectively. Denote also by vjf 1 the set of all elements in V n of 
degree < i. These subspaces define a filtration on V n : 

{0} = vt n ~ 2) c vy n -^ c • • • c vyv = v n . 

Let E be a one-dimensional vector bundle over a one-dimensional 
manifold M with a local coordinate x. Denote by J n (E) the n-th 
order jet bundle of E, which is a (n + 1) -dimensional vector bundle 
over M. Then any (n + l)-th order linear homogeneous ODE can be 
considered as a connection on J n (E) such that all its solutions, being 
lifted to J n (E), are horizontal. Let J r {J n {E)) be the frame bundle of 
J n (E). Since dim V n — n + 1, we can identify jF(J ra (i5)) as a set of all 
isomorphisms <f) x : V n — * J™(E). This turns jF(J n (i5)) into a principle 
GL(K)-bundle over M. 

Denote by uo the corresponding connection form on J r {J n {E)). In 
brief, the main idea of the Se-ashi works is that Wilczynski invariants 
can be interpreted in terms of a natural reduction of the GL(Ki)-bundle 
jF(J n (E)) to some G-subbundle P characterized by the following con- 
ditions: 

(1) G is the image of the lower-triangular matrices in GL(2,M) 
under the representation p n : GL(2,M.) — > GL(V n ); 

(2) oo\p takes values in the subspace (X, H, Z, Y, Y 2 , . . . , Y n ). 

The form uj\ p can be decomposed into the sum of ou g i with values 
in g[(2, M) C Ql(V n ) and YTi=2 UJ '^ % • Then lu s1 defines a flat projective 
structure on the manifold M, while the forms u>i (or, more precisely, 
their values on the vector field J^) coincide up to the constant with 
Wilczynski invariants 9 i+ i. 

Let us describe this reduction in more detail. Denote by G (0) the 
subgroup of GL(V n ) consisting of all elements preserving the filtration 
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Vn on V n introduced above. These are exactly all elements of GL(V n ) 
represented by lower-triangular matrices in the basis {E , . . . , E n }. For 
k > 1 denote by #ft k \V n ) the following subalgebra in gl(V n ): 

t«(K) = {0 € i(K) I Wtf>) c itf**)}. 

Let GL^{V n ) be the corresponding unipotent subgroup in GL(V n ). 
Define the subgroups G^ C G (0) as the products GL®(2, R)GL^(V n ) 
for each k > 0, where GL (0) (2,IR) is the intersection of G (0) with 
p n (GL(2, R)). Denote also by VF the subspace in gl{V n ) spanned by 
the endomorphisms Y 2 , . . . ,Y n (here Y C jj[(2,R) is identified with 
the corresponding element in Ql(V n )). 

The reduction P C JF( J n (E)) is constructed via series of reductions 
Pfe + i C Pfe, where P k is a principal G^^-bundle characterized by the 
following conditions: 

a) P consists of all frames <f) x : V n — > J£(E), which map the filtration 
of V n into the filtration on each fiber J™(E); 

b) for k > 1 the form Uk = oj\p k takes values in the subspace W k = 

At the end of this procedure we arrive at the principal bundle P = P n +\ 
with the structure group G = G n+ i and the 1-from u> = 0J n +i with 
values in W + qI(2,R). 

In fact, the second condition can be considered as a definition of Pk 
for k > 1. Indeed, let (x,y ,yi, . . . ,y n ) be a local coordinate system 
on J n (E), such that the n-jet of the section y(x) of J{E) is given by 
= y( l \x). Then the connection on J n (E) corresponding to the 
linear homogeneous equation (1) is defined as an annihilator of the 
forms: 

0i = dyi - y i+1 dx, 

n = dy n + {^2pi{x)y}j dx. 

If s: x i — ^ (y (x), . . . , y n ( x )) is any section of J n (E), then the covariant 
derivative of s along the vector field has the form: 

n 

V^s = (y' (x) -yi(x),y[(x) - y 2 (x), . . . , y' n (x) + ^p^y^x)). 

Let i — 0, . . . , n, be the standard sections defined by = 5ij. 
Then s = {so, ■ ■ ■ , s n } is a local section of the frame bundle Pq. Let 
uo be the connection form. The pull-back s*uo can be written in this 
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coordinate system as: 



S LU = 



( 







\p (x) pi{x) p 2 (x) 






-1 

Pn(x)J 



dx 



Then Se-ashi reduction theorem says that there exists such gauge trans- 
formation C: M -> G<® that 

n 

C- 1 (s*lo)C + C~ l dC = ( - X + aH + (3Z + 7 F + Oi+iY*) dx, 

and the set of such transformations forms a principal G-bundle. As 
shown in [15], the coefficients 6 i+ i coincide with classical Wilczynski 
invariants 6> i+1 defined by (3) up to the constant and some polynomial 
expression of invariants of lower weight: 

9 i+ i = c9 i+ i + P i+ i(03, . . . , 

where Pj+i is some fixed polynomial without free term. In particu- 
lar, Theorem 1 remains true if we substitute invariants 0{ with their 
modified versions Q{. 

We shall not repeat the computations from Se-ashi work [15], just 
mentioning that it is based on the following simple technical fact. 
Namely, for any k > 1 consider the subspace gi k (V n ) of all operators of 
degree k and the mapping ad k (X): gl k (V n ) — > 0[ fc „ 1 (l / „), A i— > [X, A]. 
Then we have the decomposition gl k (V n ) = (Y k ) © Imad^+iX for all 
k > 1, which allows to carry effectively the reduction from to Pk+i- 
On each step such reduction involves only the operation of solving linear 
equations with constant coefficients and differentiation. In particular, 
this proves, that Wilczynski invariants are polynomial in terms of the 
coefficients po(x), . . . ,p n (x) of the initial equation (1). 

In general, we can not make further reductions without assumption 
that some of these invariants do not vanish. And if all these invariants 
vanish, then our equation is equivalent to the trivial equation = 
0. 



4. Generalization of Wilczynski invariants to non-linear 

ODEs 

Let S be now an arbitrary non-linear ODE of order (n + 1). As 
above, we identify it with a submanifold of codimension 1 in the jet 
space J n+1 . Then in local coordinates (x, yo, ■ ■ ■ , y n +i) the equation 
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£ is given by y n+1 = f(x, y , . . . , y n ), and the functions (x, y , . . . , y n ) 
form a local coordinate system on £ . 

Let us recall that there is a canonical contact C k distribution defined 
on each jet space J h (R 2 ), which is generated by tangent lines to all lifts 
of curves from the plane. In local coordinates it is generated by two 
vector fields: 

/ d d d d 

C = (^— , ^ + Vi^ + ■■■y k - 



dy k dx dy dy k -i 

We also have natural projections n k ^: J k {R 2 ) — > J'(R 2 ) for all I < k. 

The contact distribution C n+1 defines a line bundle on the equation £ 
as follows E p = T p £ n C r p l+1 for all p G £. Integral curves of this bundle 
are precisely lifts of the solutions of the given ODE. In local coordinates 
this line bundle is generated by the vector field: 

„ d d d „ d 



dx dy dy n _i dy n 

which defines also the total derivative operator. 

We are interested in the frames on the normal bundle to E, that is 
the bundle N(£) = T£/E. We have a natural filtration of this bundle 
defined by means of the projections n n+lti : £ — > J l (R 2 ) for i < n. 
Namely, we define iVj as the intersection of vr* +1 { C l with T£ modulo 
the line bundle E for all i — 0, . . . , n. Then it is easy to see that the 
sequence 

N = N D JVi D . . . N n D 0, 

is strictly decreasing and dim A^j — n + 1 — i. In local coordinates we 
have 

*-<£--£> + * 

The vector field D defines a first order operator on N, which is com- 
patible with this filtration, i.e. D(Ni) C iVj_i for alH = 1, . . . , n. 

As in the case of linear equations, we can define the frame bundle, 
consisting of all maps (f) p : V n — > N p from the GL(2 : i?)-module V n into 
N p , which preserve the nitrations. This defines a G^-bundle P over 
£. 

To proceed further with similar reductions, we need to have an analog 
of the connection form. In general, we don't have any natural connec- 
tion on the normal bundle N(£). However, we can define the covariant 
derivative along all vectors lying in the line bundle E generated by D: 

V D (X + E) = [D, X] + E, for any X + E G N(£). 
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This gives us a so-called partial connection with the connection form 
uj: E — > gl(V n ). This connection form is sufficient for our purposes. As 
in the case of linear equations, we can form the similar set of reductions 
Pi of the principle bundle Pq. Each such bundle Pj, being restricted to 
any solution L, coincides with the principle bundle constructed from 
the linearization of £ along this solution. As a result, we get the set of 
well-defined generalized Wilczynski invariants Wi, where each Wj is a 
global section of the bundle S l E* . 

Since the bundle N(£) and the connection form uj are defined to- 
tally in terms of the contact geometry of the jet space J n+l and the 
reduction procedure does not depend on any external data, we see that 
generalized Wilczynski invariants are contact invariants of the original 
equation. It can be formulated rigorously in the following way: 

Theorem 3. Assume that n > 2. Let <p: J n+1 — > J n+1 be a local con- 
tact transformation establishing the local equivalence of two equations 
£ and £. Let E (E) be the line bundle on £ (resp. £) defining the so- 
lution foliation. Let Wi G T(S l E*) (resp. Wi G T(S l E ) be generalized 
Wilczynski invariants of £ (resp. £). Then we have <fi*(E) = E and 
<P*Wi = Wi for alli = 3,...,n + l. 

This theorem was first proven in the work [5] via direct computa- 
tion of the canonical Cartan connection associated with any ordinary 
differential equation of order > 4 (see [17, 4] for more details). 

5. Equations with vanishing Wilczynski invariants 

5.1. Structure of the solution space. Suppose now that all gener- 
alized Wilczynski invariants vanish identically for a given equation £. 
For linear equations this would mean that the given equation is triv- 
ializable, and the associated curve L £ is an open part of the rational 
normal curve in IRP n . 

Even though the curve L £ is defined up to contact transformations, 
we can define geometric structures on the solution space itself, which 
will be independent of the choice of fundamental solutions. Indeed, let 
V(£) be the solution space of the linear equation (1). Then at each 
point x on the line we can define a subspace V xo (£) in V(£) consisting 
of all solutions vanishing at x . Thus, assigning the line Vj^ in the 
dual space V(£)* to each point xq we get a well-defined curve in the 
projective space PVg . 

Using the duality principle, we can also construct the curve in the 
projective space PVg. To do this, we can consider the line l Xo (£) in 
V{£) consisting of all solutions vanishing at x with their derivatives 
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up to order n — 1. The mapping x \— > l Xo (£) defines a non-degenerate 
curve in the projective space PV(£). 

Suppose now, that all Wilczynski invariants of the linear equation £ 
vanish identically and consider the symmetry algebra g C sl(V(£)) of 
the constructed curve. Fixing any fundamental set of solutions and ap- 
plying Theorem 1, we can show that this Lie algebra is isomorphic to 
sl(2, M) and acts irreducibly on V{£). Let G be the corresponding sub- 
group in SL(V(£) and let G be its product with the central subgroup 
in GL(V(£)). Then G is naturally isomorphic to GL(2,WL), and the 
G-module V{£) is equivalent to the GL{2, M)-module S n {R 2 ). Thus, 
any equation with vanishing Wilczynski invariants defines a CL(2,R)- 
structure on its solution space. 

For a non-linear equation the vanishing of generalized Wilczynski 
invariants is not sufficient for being trivializable. The paper [5] de- 
scribes the extra set of invariants that should vanish to guarantee the 
trivializability of the equation: 

Theorem 4 ([5]). The equation (4) is contact equivalent to the trivial 
equation = if and only if its generalized Wilczynski invariants 

vanish identically and in addition: 

forn = 2: f 222 2 = 0; 

for n = 3: / 333 = 6/ 2 33 + /f 3 = 0; 

for n = A: f u = 6/ 234 - 4/ 333 - 3/f 4 = 0; 

forn = 5: / 55 = / 45 = 0; 

for n>6: f n>n = f n ,n-l = /n-l,n-l = 0. 

Yet, even if the equation is not trivializable, its linearization at each 
solution is trivializable, and we can construct a family of associated ra- 
tional normal curves. It will define a GL(2, IRQ-structure on the solution 
space S, if it is a Hausdorff manifold. 

Let £ be an arbitrary non-linear ODE. Suppose that its solution set 
S is Hausdorff (i.e., there exist a quotient of £ by the foliation formed 
by all solutions). Let yo(x) be any solution of the non-linear ODE. This 
is just a point in the manifold S. The tangent space T yo S to S at yo(x) 
can be naturally identified with a solution space of the linearization 
of £ along yo(x). By the above, we can naturally construct a curve / in 
the projectivization of T yo S, which means that we have a well-defined 
two-dimensional cone C yo in each tangent space T yo S. 

Suppose now that all generalized Wilczynski invariants of £ vanish 
identically. Then linearizations of £ along all solutions are trivializable, 
and all cones C yo C T yo S are locally equivalent to the cone in R ra+1 
corresponding to the normal curve in RP n . 

Thus, we arrive at the following result. 
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Theorem 5. Let £ be an arbitrary (non-linear) ODE with vanishing 
generalized Wilczynski invariants. Suppose that its solution space S 
is Hausdorff and, hence, is a smooth manifold. Then there exists a 
natural irreducible G L (2, R)- structure on S. 

This structure was constructed in [6] in a slightly different way and 
is called a paraconformal structure. 

5.2. Examples from twistor theory. In general, explicit examples 
of equations with vanishing Wilczynski are very difficult to construct. 
However, as suggested by N. Hitchin [11], there is a large class of ex- 
amples coming from twistor theory. 

The whole theory above is also valid in complex analytic category. 
Consider an arbitrary complex surface S with a rational curve L on 
it. Suppose that the normal bundle of L has Grothendieck type 0(n). 
Then Kodaira theory states that this rational curve L is included into 
a complete (n + l)-parameter family £ = {L a } of all deformations of 
L. 

The family C uniquely defines an ordinary differential equation of 
order (n + 1) such that all these curves are its solutions. In more 
detail, we have to consider the lifts of the curves L a to the jet space 
J n+1 (S), where they will form a submanifold of codimension 1. This 
defines an ordinary differential equation £. 

Theorem 6. Let £ be an ordinary differential equation defining the 
complete set of deformations of a rational curve L on a complex surface 
S. Then all generalized Wilczynski invariant of £ vanish identically. 

Proof. Each generalized Wilczynski invariant Wj, being restricted to 
the lift of any solution L a , defines a global section of the line bundle 
S l T*L a , which has Grothendieck type 0{— 2i). Since the only global 
section of the line bundle of this type is zero, we see that all generalized 
Wilczynski invariants will vanish identically. □ 

Consider a number of explicit examples. 

Example 1. Let S = CP 2 and let L be a quadric in S. Then the 
normal bundle of L has type 0(4), and the complete family of defor- 
mations is 5-dimensional. Clearly, all quadrics in CP 2 belong to this 
family and depend on exactly 5 parameters. 

So, we see that the family C in this case is a family of all quadrics on 
the complex projective plane. The differential equation of all quadrics 
is well-known and has the form: 



9(/) 2 y( 5 )-45 2/ V , l/ (4) +40( 2/ w ) 3 = 0. 
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By above, we have W3 = W4 = W5 = for this equation. However, 
this equation is not trivializable, since, for example, its symmetry al- 
gebra is only 8-dimensional, while the trivial equation = has a 
9-dimensional symmetry algebra. 

Example 2. Let S be a rational surface Sk = P(O(0) + 0(k)) viewed 
as a projective bundle over CP 1 . Then all sections of this bundle with 
a fixed intersection number I with fibers form a complete family of 
deformations depending on k + I + 1 parameter. In the appropriate 
coordinate system on Sk they can be explicitly written as: 

a + aix H h a k x h 

V{X) ~ b + b lX + --- + kx l ' 

Assume that k >l. (Otherwise we can substitute y with 1/y.) This 
set of rational curves defines the following differential equation: 

= 0, where Zi = - — r~y®. 

For example, in the simplest case I = we get just the trivial equa- 
tion = 0, while for I — 1 we get the following equation: 

(k + l)y^y ik+ V - (k + 2)(y^) 2 = 0. 

For k > 2 the symmetry algebra of this equation is k + 2-dimensional, 
and, hence, it is not trivializable. 

For I > 2 the equation F^i = is never trivializable, since, for 
example, being solved with respect to the highest derivative y( fc+z+1 ) it 
is not linear in term y( k+l \ 

Since all solutions of the equation F k j = are by construction ratio- 
nal curves, all its generalized Wilczynski invariants vanish identically. 

Example 3. Consider the standard symplectic form a on the 4- 
dimensional vector space C 4 . Then it defines the contact structure 
on CP 3 invariant with respect to the induced action of SP(4, C). It is 
easy to see that the set of all rational normal curves in CP 3 , which are 
at the same time integral curves of the contact distribution, depends 
on 7 parameters and forms the complete family of deformations of any 
of such curves. 

The equation describing this set is given by: 

10(y'")V 7) - 70(y'")V 4 y 6) - 4%'")% (5) ) 2 

+ 28(W 4) )y 5) - 175(y (4) ) 4 = 0. 



Zk-l+l Zk-l+2 ■ ■ ■ Zk+1 
Zk-l+2 Zk-l+3 ■ ■ ■ Zk+2 

z k+l z k+2 ■ ■ ■ z k+l+l 
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Again, since all solutions of this equation are rational curves, all its 
generalized Wilczynski invariants vanish identically. Yet, this equation 
is not trivializable, since its symmetry algebra is 10-dimensional (in 
fact, it coincides with sp(4, C)), while the trivial equation = has 
11-dimensional symmetry algebra. 

Note that equations from the examples above have the common 
property that they have the symmetry algebra of sub maximal dimen- 
sion [13]. 



6. Wilczynski invariants for systems of ODEs 

Since all results of this paper are based on two ideas, namely the 
notion of linearization and the invariants and structures associated with 
linear equations, they are directly generalized to systems of ordinary 
differential equations and even to equations of finite type. 

The generalization of Wilczynski invariants to the systems of ODEs 
was obtained by Se-ashi [16]. Consider an arbitrary system of linear 
ordinary differential equations: 

y (n+l) + Pn ( x )y(n) + ... + P ( X )y( X ) = , 

where y(x) is an R m - valued vector function. The canonical Laguerre- 
Forsyth form of these equations is defined by conditions P n = and 
trP„_! = 0. Then, as in the case of a singe ODE, the following expres- 
sions: 
(6) 

k ~2-f j {k - - 1)! F n-k+p K-2,...,n + L. 

are the End(R m )-valued relative invariants, where each invariant 8{ has 
weight %. Note that unlike the case of a single ODE, the first non-trivial 
Wilczynski invariant has weight 2. 

We can also consider a characteristic curve associated with any 
linear system, which will take values in the Grassmann manifold 
Gr m (R( n+1 ) m ). Se-ashi also proved that the system is trivializable if 
and only if all these invariants vanish identically, or, equivalently, if 
the characteristic curve is an open part of the curve defined by a triv- 
ial equation. The symmetry group of this curve is isomorphic to the 
direct product of SX(2,R) and GL(m, R) with an action on R(" +1 ) m 
equivalent to the natural action on S , ^ n+1 - ) (R 2 ) <S> R m . 

All that leads us immediately to the following generalization of The- 
orems 2 and 5. 
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Theorem 7. Let £ be an arbitrary (non-linear) system of m ordinary 
differential equations of order (n + 1). Then there exist generalized 
Wilczynski invariants W2, ■ ■ ■ , W n+ i that, being restricted to each so- 
lution, coincide with Wilczynski invariants of the linearization along 
this solution. 

Suppose that all these invariants vanish identically and the solution 
space S is a smooth manifold. Then there exists a natural irreducible 
SL(2,R) x GL(m, R) structure onS. 

Remark 1 . Let us note that in the case of systems of ODEs the invariant 
Wi is a section the vector bundle S l E* <g> End(V^), where V n is a sub- 
bundle of the normal bundle N(S) = TS/E defined as < +ljn C n n T£ 
modulo E. In case of a single ODE the bundle V n is one-dimensional, 
and End(V^) can be identified with the trivial bundle over £. 

Example 4. Consider the simplest non-trivial case of second order 
systems of ODEs. The linear homogeneous system can be written as: 

(7) y" = A(x)y> + EJ(x)y, 

where y(x) is the unknown R m -valued function and A(x),B(x) G 
End(M m ). The set of all transformations preserving the class of such 
systems has the form (x,y) 1— > (X(x) , /i(x)y) , where X(x) is a local line 
reparametrization and n(x) G GL(m,R). 

We can bring the equation (7) into the semi-canonical form with 
vanishing coefficient A(x) by means of a certain gauge transforma- 
tion (x,y) 1 — (x,fji(x)y). Then, applying a reparametrization (x,y) 1— > 
(X(x),y) and an appropriate gauge transformation to bring the equa- 
tion back to the semi-canonical form, we can make trace of B(x) vanish. 

The traceless part of the coefficient B(x) in the semicanonical form 
gives us the Wilczynski invariant 2 for linear systems of second order. 
Explicitly, it can be written as: 
(8) 

© 2 = - l/mtr$(x), where <b{x) = B(x) A'(x) + -A(x) 2 . 

According to Se-ashi [16] , this is the only Wilczynski invariant available 
for systems of second order, and the equation (7) is trivializable if and 
only if the invariant ©2 vanishes identically. 

Consider now an arbitrary non-linear system of second order: 

y" = fi(x,yj,y' k ), i = i,...,m. 

Then its generalized Wilczynski invariant can be obtained from (8) by 



substituting A(x) with the matrix ( ^4- ) , B(x) with the matrix 
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and the usual derivative d/dx with operator of total derivative: 



Denote by Wi the invariant we get in this way. It coincides (up to the 
constant multiplier) with the invariant constructed by M. Fels [7]. 

If the invariant Wi vanishes identically, then we get the GL(2, R) x 
GL(m, R) structure on the solution space. This structure is also called 
Segre structure and was first constructed for solution space of sec- 
ond order ODEs with vanishing generalized Wilczynski invariant by 
D. Grossman [9]. 
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